High-dimensional fractionalization and spinon deconfinement in pyrochlore antiferromagnets 
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The ground states of Klein type spin models on the pyrochlore and checkerboard lattice are spanned by the 
set of singlet dimer coverings, and thus possess an extensive ground-state degeneracy. Among the many exotic 
consequences is the presence of deconfined fractional excitations (spinons) which propagate through the entire 
system. While a realistic electronic model on the pyrochlore lattice is close to the Klein point, this point is in 
fact inherently unstable because any perturbation e restores spinon confinement at T = 0. We demonstrate that 
deconfinement is recovered in the finite-temperature region e <C T <C J, where the deconfined phase can be 
characterized as a dilute Coulomb gas of thermally excited spinons. We investigate the zero-temperature phase 
diagram away from the Klein point by means of a variational approach based on the singlet dimer coverings of 
the pyrochlore lattices and taking into account their non-orthogonality. We find that in these systems, nearest 
neighbor exchange interactions do not lead to Rokhsar-Kivelson type processes. 

PACS numbers: 75.10.-b, 75.50.Ee, 75.40.Cx, 73.43.Nq 



I. INTRODUCTION 

Discovering new phases of matter is a primary objective 
of physics. The fractionaUzed spin liquid in two spatial di- 
mensions {d = 2) has provided a popular candidate frame- 
work for models describing the exotic properties observed in 
many strongly correlated electronic materials, including frus- 
trated quantum magnets and high-Tc superconductors lU. 
The spin-liquid state is characterized by spinon excitations 
carrying unit charge under a compact U( l) gauge field. How- 
ever, Polyakov has argued |2f| that a pure compact U(l) gauge 
theory is always confining at zero temperature for d ~ 2, con- 
finement between test particles with opposite charges being 
produced by the proliferation of instanton tunneling events. 

By contrast, the case for confinement by instantonproUf- 
eration in spin systems is rather more involved iH Ql, and 
the situation at a critical point offers additional possibilities. 
These considerations underlie a recent discussion of the de- 
confined quantum critical point jst] as the foundation for the 
emergence of spin-liquid phases with fractional excitations. 
However, while this set of elegant ideas appears most plausi- 
ble, it is not yet clear from any analytical or numerical calcula- 
tions that there exists a microscopic Hamiltonian, and particu- 
larly one with only simple electronic interactions, exemplify- 
ing this class of physical phenomena. Quasiparticle fraction- 
alization has been demonstrated rigorously in one dimension 
(Id) and in the quantum Hall effect. For frustrated quantum 
spin systems, recent progress was made in the context of a 
two-dimensional model with high degeneracy in the ground- 
state manifold which exhibits (partial) dimensional reduction 
and consequent spinon deconfinement We will expand 
upon and systematize these ideas to provide a further demon- 
stration of a class of physical systems displaying fractional 
excitations. 

To this end, we begin with the natural question of whether 
a dilute gas of deconfined spinons can be stabilized in dimen- 
sions d > 1 in a finite region of the phase diagram. This 
latter requirement is crucial for observing deconfined spinon 



excitations in real systems. We will provide an affirmative 
answer by considering a highly frustrated 5 = 1/2 model 
on the d ~ 2 (checkerboard) and d = 3 pyrochlore lattices. 
We demonstrate that for T = the spinon excitations are de- 
confined at the Klein critical point, whose nature we will ex- 
plain in detail. As expected from above, this zero-temperature 
deconfined phase is unstable under any realistic perturbation 
characterized by an energy scale e. We argue that in the regime 
e ^ J, where J is the characteristic energy scale of the model, 
a dilute Coulomb gas of spinons is stabilized in the tempera- 
ture range e ^ T ^ J. This result raises the possibility of 
observing spinons in real higher-dimensional systems. 

The checkerboard and pyrochlore lattices (Fig. [TJ which 
form the focus of our analysis have long been known to pos- 
sess a geometry which is amenable to very large ground- 
state degeneracies. The pyrochlore lattice is composed of 
corner-sharing tetrahedra, arranged in a cubic (fee) struc- 
ture, such that each vertex is common to two 4-site units 
(see panel (b) of Fig.[T]i. This geometry is rather common in 
transition-metal and rare-earth oxide systems, specifically in 
the pyrochlore (the origin of the name) and spinel structures. 
The pyrochlore systems Ht] have chemical formula A2B2O7, 
where both A and B may be taken from a wide range of metal 
ions, particularly rare-earths, and both occur on a pyrochlore 
lattice. Spinel systems have chemical formula AB2O4, 
where it is more typical for the A and B ions to be transi- 
tion metals and only the B sites form a pyrochlore structure. 
Individual compounds in these classes are usually strongly in- 
sulating and of fixed ionic valences, and the presence of non- 
magnetic ions on one of the sublattices (the A sites for the 
spinel) leads to a pyrochlore lattice of interacting spins on the 
other. It was noted at a very early stage ifiol [Till that the con- 
nectivity of such a structure is not dissimilar to the arrange- 
ment found in water ice, a point which we will both explain 
and exploit extensively below. 

The checkerboard lattice can be considered as a 2d version 
of the pyrochlore structure, in the sense that it also is com- 
posed of corner-sharing tetrahedral units, these being repre- 
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sented by the crossed plaquettes (the terminology we adopt 
henceforth) in panel (a) of Fig.[T] It is the relatively low coor- 
dination number and tetrahedral connectivity of the checker- 
board and pyrochlore lattices which permits a proliferation of 
ground states. Some of the essential considerations are pre- 
sented for classical pyrochlore antiferromagnets in Ref. lll2ll 
and for quantum pyrochlore antiferromagnets in Ref. lll3ll . 
While a significant body of recent work has led to a better 
understanding of ice-like phases in pyrochlore ferromagnets 
lfl4ll . the situation for both classical and quantum antiferro- 
magnets remains far from satisfactory, to the extent that cer- 
tain results conflict lITsll . Here we provide a further contribu- 
tion to the body of knowledge concerning quantum pyrochlore 
antiferromagnets . 

We will demonstrate further that the highly unconventional 
properties listed above arise within simple, nearest-neighbor, 
bilinear and biquadratic spin models on the 2d and 3d py- 
rochlore lattices. These models are derived directly from 
half-filled, nearest-neighbor Hubbard models, i.e. from mod- 
els with Coulomb interactions of the shortest range, and are 
thus fully realistic for electronic systems. The key features 
giving rise to exotic behavior are (i) an exponential number 
of degenerate ground states; (ii) exact critical behavior with 
an emergent divergent correlation length, accompanied by the 
absence of long-ranged order; (iii) fractionalization of excited 
spin states into spinons which propagate freely in space, im- 
plying a separation of spin and charge upon doping these sys- 
tems (we will show that the emergent fractionalized excita- 
tions have well-defined quantum numbers which may be de- 
scribed as residing at the ends of strings obeying simple rules); 
(iv) an effective dimensional reduction occurring precisely at 
the critical point; (v) topological order, in a sense which we 
will define explicitly below. All of these features lead in com- 
bination to a distinct type of critical point whose static and dy- 
namic properties we will characterize in detail, illustrating in 
addition their connection to a multitude of ideas and concepts 
emerging in the recent literature in the context of deconfine- 
ment and criticality. 

In Sec. II we introduce the type of spin model display- 
ing a Klein point, and discuss its foundation in the electronic 
Hamiltonian of real, correlated insulating materials. In Sec. Ill 
we begin our analysis of the model at the Klein critical point 
and at zero temperature, focusing on the characteristic proper- 
ties here and illustrating their origin in the extensive ground- 
state degeneracy. In Sec. IV we remain at the Klein point but 
consider finite temperatures, to derive the effective entropic 
interaction of spinon excitations and thus demonstrate their 
thermally driven deconfinement. In Sees. V and VI we restore 
some of the perturbing effects present in the realistic Hamil- 
tonian of Sec. II. Working at zero temperature, in Sec. V we 
introduce and apply a specifically constructed variational pro- 
cedure based on the (non-orthogonal) dimer coverings of the 
ground-state manifold, to demonstrate the selection of spe- 
cific valence-bond orderings as the ground states of the per- 
turbed system. In Sec. VI we extend our considerations to fi- 
nite temperatures and present one example of an exactly crit- 
ical regime. These calculations indicate explicitly both the 
confinement of spinons at low temperatures and their liber- 
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FIG. 1 : (Color online.) (a) Checkerboard and (b) pyrochlore lattices. 



ation at higher temperatures, thus reinforcing the results of 
Sees. Ill and IV. Section VII contains a discussion and con- 
clusion. 



II. MODEL 

The spin model upon which we focus is a = 1/2 Hamil- 
tonian on the checkerboard and pyrochlore lattices. 



-K^{H^^H^i+HSHf, + H^,Hfi), (1) 



h = jJ2h^- 

where (ij) denotes all pairs of sites in the tetrahedron a = 
ijkl and H^j = ■ S" is a scalar product of the two spins. 
On the checkerboard lattice a denotes each crossed plaque- 
tte, by which is meant those in Fig.lTta) with cross-coupling 
interactions, and the first term is equivalent to both nearest- 
and next-neighbor Heisenberg interactions of strength J. This 
model was considered on the square lattice (all plaquettes 
crossed) by two of the present authors I?!], and we will con- 
trast the two d = 2 systems below. The Heisenberg interaction 
on the tetrahedral units (we will also refer to the crossed pla- 
quettes of the checkerboard lattice as tetrahedra) of the d = 2 
and d = 3 pyrochlore lattices is strongly frustrated, and the 
symmetric four-spin interaction pushes the model to maximal 
frustration at Kc = 4 J/5 ifsoll. 

A straightforward but important observation is that H may 
be recast in the form 



with the net spin of the tetrahedral unit M, 

= ^ Si, 



(2) 



(3) 



Ji = J - 7K/A and J2 = K/2 tM- This Hamiltonian O 
can be obtained from the Hubbard model on the pyrochlore 
lattice, 

fl"Hubb = X! '^L^jcr + U ^n^in^, (4) 
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where riia = c\^Ci„ is the number operator for electrons of 
spin a. This is perhaps the simplest possible description of in- 
teracting electrons with only nearest-neighbor hopping t and 
a local Coulomb repulsion U . Systems which are half-filled 
and have strong on-site interactions are localized insulators 
where, with the exception of virtual processes, the kinetic en- 
ergy gain from hopping is sacrificed to avoid the Coulomb 
cost of double occupancy. The effective Hamiltonian obtained 
from the virtual hopping processes at lowest (second) order 
in the small parameter t/U is a nearest-neighbor Heisen- 
berg model of the type contained in the first term of Eq. 



with Jl = At^ /U and electronic spin 5f 

(a = X, y, z). 

A fourth-order strong-coupling expansion of Eq. ^ re- 
turns interaction terms bilinear in spins on the same and on 
different tetrahedra and biquadratic in spins on the same tetra- 
hedron 1I32I1 . As shown in detailed studies performed for pla- 
quettes of the square lattice in the context of planar cuprates 
II33I1 . and verified by fits to magnon dispersion relations for 
the same systems 13411 . the latter interactions are stronger than 
the former by approximately one order of magnitude. Specif- 
ically, a far larger number of intra-tetrahedron processes, 
whose sum has the symmetric form contained in the sec- 
ond term of Eq. ([T]i, contributes to the effective fourth-order 
Hamiltonian than do processes of the next-neighbor Heisen- 
berg type. Thus by retaining only the intra-tetrahedron pro- 
cesses to fourth order in t one obtains the spin Hamiltonian of 
Eq. © with iH 
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(5) 



The general, SU(2)-invariant spin Hamiltonian arising from 
a realistic electronic model may thus be considered as an 
intra-plaquette interacting system of exactly the form given 
in Eq. (|2]l, with only very weak perturbations from inter- 
tetrahedron terms. To quantify this statement, in the half- 
filled system Eqs. (|2|i and ^ are related by 



Hubb 



Si • Sj, 



(6) 



where i?Hubb is the effective spin Hamiltonian obtained from 
a fourth-order strong-coupling expansion of i^Hubb- The sum 
is over pairs of sites {{i, j)) separated by two bonds and which 
do not belong to the same tetrahedron. The exchange coupling 
J3 is given by 



It is evident from Eqs. (|5]i and (|7]i that for small t/U the 
inter-plaquette corrections ( J3) originating from the Hubbard 
Hamiltonian are indeed one order of magnitude smaller than 
the intra-plaquette terms of H A particle-hole transfor- 
mation leaves the spectrum of i^Hubb invariant at half-filling. 



Since the net effect of this transformation is to change the sign 
of the hopping, iJnubb cannot contain any term of odd order 
in t. 

In Sec. Ill we wiU analyze the properties of Eq. (|2|i at 
the Klein point, Ji = -J2 H or K ^ Kc ^ 4 J/5, 
where the semi-positive definite Hamiltonian of Eq. ^ at- 
tains its minimal value (when 5*^1 = or 1 on every tetrahe- 
dron). The Klein point is quasi-exactly solvable. [Its ground 
states may be determined exactly.] Following the observations 
in the preceding paragraph, in Sees. V and VI we will ex- 
pand our treatment to include perturbations around the Klein 
point, which are of two types: the stronger, represented by 
nearest-neighbor Heisenberg terms, correspond to deviations 
of Eq. ([T]) from the point K = Kc', the weaker, represented 
by next-neighbor Heisenberg terms (and, for the checker- 
board, cyclic terms on the uncrossed plaquettes), correspond 
to fourth-order inter-plaquette interactions. Specific S* = 1/2 
materials can be expected at least in the 3d pyrochlore struc- 
ture, such as that formed by the magnetic ions at both A'"+ 
and B"+ {m+n = 7) sites in the stereotypical pyrochlore sys- 
tems A2B2O7, and on the B sites in A"'+B2 +O4 (m+2?7, = 8) 
spinels. However, despite the considerable choice of magnetic 
and non-magnetic ionic species afforded by these two struc- 
tural classes, to date we are not aware of the successful syn- 
thesis of any such compounds which maintain the full cubic 
lattice symmetry. 



III. EXACT GROUND STATES FOR THE 5 = 1/2 SYSTEM 
AT THE KLEIN POINT 



Equation Q suggests an even simpler way of expressing 
the Hamiltonian, 



(8) 



with a quartic polynomial in the total spin of the four 
sites forming a given tetrahedral unit of the pyrochlore lattice, 
hm = \JiS'^ + \ J2S^. For clarity we provide a brief re- 
view of the possible two- and four-particle states of = 1/2 
spins. For a system of four coupled spins, the total spin of 
any tetrahedron obeys < 5*^ < 2. The 16-state spin space 
is decomposed into two singlet states iS^ = 0), three triplets 
(5^ = 1), and one quintet (S*^ = 2), 



1 1 
- ® - 

2 2 



0©0® 1®1©1^ 



(9) 



In this explicit decomposition the right-hand side labels the 
disjoint net spin (5^) sectors while the left-hand side encodes 
the 16-dimensional space spanned by the direct product of the 
four S* = 1/2 particles. The sum of any two nearest-neighbor 
spins on the lattice [Si + Sj)"^ = Spairl^pair + 1) takes only 
the two values 5pair = 0, 1. If any two spins within a given 
tetrahedron are in a singlet state (5pair = 0), then the total 
spin of this tetrahedron cannot exceed unity, S*^ < 1. For 
the parameter choice K = Kc ( J2 = — Ji), Jf (Eji can be 
rewritten, up to irrelevant constants, as a Klein Hamiltonian 
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FIG. 2: (Color online.) Highly regular ground states on the checker- 
board and pyrochlore lattices. The ovals denote singlet dimer states. 
The arrows denote the representations of these dimer states within 
the six-vertex model (see text). On each plaquette (tetrahedron) the 
dimer connects the bases of the two incoming arrows. 
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(10) 



where V is the projection operator onto the subspace of net 
spin S^ = 2. 

Next we note, following Ref. ll39ll . that because the number 
of lattice sites is double the number of tetrahedra, TVs = 2A^t, 
the set of ground states with one singlet per tetrahedron can 
be mapped onto the set of hard-core dimer coverings of the 
lattice. Thus any dimer covering of the lattice with precisely 
one dimer per tetrahedron is a ground state of Hk- Further, 
any state \ tp) which is a superposition of dimer coverings each 
of which has one dimer per tetrahedron. 



p 



ijeP 

(11) 

(the dimer coverings are labeled by P), is also a ground state 
of Hk- One of the simplest ground states is afforded by the 
dimer coverings depicted in Fig.|2] While it is clear that any 
state of the form of Eq. (fTTT i is a ground state, a proof that 
all ground states are of this form is far less obvious. It may 
nevertheless be shown rigorously [[ssll that for Klein models 
on the pyrochlore lattices all ground states are indeed of the 
type specified in Eq. ( fTTT i. 

As a first qualitative consequence of this result, it was noted 
in Ref. ff] for the analogous point in the square-lattice Hamil- 
tonian that satisfying this constraint leads to a significant de- 
generacy because dimer singlets may be rearranged along di- 
agonal lines of the lattice with no energy cost. The ground- 
state entropy scales with the perimeter of the system, and pairs 
of single-spin (spinon) excitations on the diagonal lines have 
no binding energy for any separation. Thus the effective di- 
mensional reduction from 2d to Id leads to the presence of 
deconfined spinons whose propagation is essentially unidi- 
mensional. On the pyrochlore lattice the open structure of 



crossed plaquettes creates a less constrained system, and the 
number of dimer configurations in the ground-state manifold 
is strongly enhanced, the still more massive degeneracy im- 
plying similarly exotic physics in this case. We proceed quan- 
titatively by computing the ground-state entropy in Sec. IIIA. 

In addition to the simple coverings shown in Fig. |2l a far 
richer variety of states exists. The mapping onto the spin-ice 
problem mentioned briefly above provides a useful classifica- 
tion of these ground states. We stress that the states with one 
dimer per tetrahedral unit defined by Eq. (fTTT ) map exactly to 
the spin-ice problem ifTil [l6ll where, motivated by the struc- 
ture of H atoms in solid water ifTolfTlll . two of the sites of any 
elementary unit (a tetrahedron) are associated with an ingoing 
arrow pointing towards the center of the tetrahedron and two 
sites lie on arrows pointing outwards. We label the two sites 
belonging to a singlet dimer | Sij ) in a given tetrahedron Kl by 
two incoming arrows from sites i and j to the center of the unit 

= ijkl. In this fashion it is clear that the system is mapped 
to a set of continuous directed lines such that each tetrahedral 
unit has exactly two incoming and two outgoing arrows rel- 
ative to its center. An example of this mapping is illustrated 
in Fig. |2] A longer but ultimately equivalent version of this 
mapping, enabling a study of correlations and degeneracies, 
was reported in Ref. js^]. The mapping is one-to-one: any 
spin-ice configuration determines a unique singlet-covering 
state with one dimer per tetrahedral unit and vice versa. The 
conservation of incoming and outgoing arrows in the spin-ice 
representation is the foundation for the divergence-free con- 
dition we will employ below. Throughout this work we fo- 
cus on systems with periodic boundary conditions (PBCs). A 
system with open boundary conditions (OBCs) would possess 
a number of additional states not captured by the six-vertex 
mapping, the entropy of which scales with its surface. 



A. Ground-State Degeneracy and Entropy 

Because the number of pyrochlore spin-ice states is 
bounded from below i40i1 by the Pauling limit illil . the num- 
ber of ground states (or, more precisely, of non-orthogonal 
dimer coverings) is given by 



Ng > (3/2)^/^ 



(12) 



with N the number of vertices of the pyrochlore lattice and 
Ng the number of ground-state singlet dimer coverings. We 
emphasize that this degeneracy is exponential in the system 
volume, yielding an extensive entropy with a simple, ana- 
lytical expression for the lower bound, the Pauling entropy 
Sp = ^NkB In |. The "exact" entropy of the ice problem is 
known from series expansion methods ll4lll . which yield the 
numerical result 

Sg ^ (0.20501 . . . )NkB > Sp = (0.20273 . . . )NkB. (13) 

Similarly, for the checkerboard lattice the system is reduced 
to the exact spin-ice model ll42ll on the square lattice formed 
by the centers of the crossed plaquettes, leading to an exact 
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FIG. 3: (Color online.) Standard representation of six-vertex states 
in terms of lines ll42ll . Every line is composed of links whose arrows 
point to the right in the vertex representation. 



(14) 



singlet dimer covering degeneracy 11421 14311 of 



with N the number of vertices of the checkerboard lattice, 
whence Sg = jNk-Q In |. In this derivation we have assumed 
that the singlet dimer coverings, although non-orthogonal, are 
linearly independent 1441 . There are additional "defect" states 
in the ground-state manifold, in which two singlet dimers may 
occupy a single tetrahedron, or un pair ed spins may be present, 
without altering the total energy ll46ll . However, because the 
entropy of these states scales with the system surface they con- 
stitute a set of measure zero in comparison with the extensive 
set of dimer coverings, and will be neglected in what follows. 



B. Geometry of ground states 

To examine the structure of general ground states of the 
model ([T]i at the Klein point, we invoke the line representation 
of the six-vertex model |42] depicted in Fig. [3] In this con- 
vention, every line is composed of links whose arrows point to 
the right, and there is clear one-to-one mapping between the 
six-vertex and the line representations. An essential property 
of the lines is their "chiral" nature: they always represent a 
motion only to the right (Fig.O. We will show in detail below 
(Fig- 13 that this chirality is not merely an artificial feature of 
the line representation, but encodes the physical restrictions 
on the possible spinon paths in the background of the allowed 
dimer states. 

Finding the elementary process connecting a given dimer 
covering with other ground states is straightforward in the line 
representation. Fig. |4|c) shows the line representation of the 
ground state depicted in Fig. |2] The simplest process in the 
dimer basis corresponds in the line representation to flipping 
over one of the line corners, as depicted in Figs.|4|c) and (d). 
Note that there is an uncrossed plaquette (in the original lat- 
tice) associated with each such corner flip. Returning to the 
dimer representation, such a process corresponds to a cyclic 
rotation of the four dimers connected to the uncrossed pla- 
quette [Fig.Ub)]. A similar process around a hexagon whose 
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FIG. 4: (Color online.) Left: representation of six-vertex states in a 
highly regular dimer configuration (a) in terms of lines (c). Joining 
the vertex coverings according to the standard prescription (see text) 
(4^1 allows one to label all permitted ground states. Right: an ele- 
mentary local spin interaction process, leading to a new dimer con- 
figuration (b), corresponds in the line representation (d) to flipping 
one of the line corners. 



edges belong to six different tetrahedra is obtained for the py- 
rochlore lattice. For clarity of explanation we comment that 
in Figs. Uc) and (d), and later in this section, the six-vertex 
(thin lines with arrows) and line representations (thick) are 
superposed to aid in illustrating their equivalence, but either 
alone is a complete representation of the system in terms of 
independent sets of quantum numbers. 

The mapping to the line representation shows that the low- 
energy sector of a Hamiltonian at the Klein point can be 
mapped into a string model. It is clear that the total num- 
ber of lines is a conserved quantity under any local physical 
process, where by "local" we mean that the process changes 
only a finite number of (spin) degrees of freedom. The dimer 
coverings may therefore be classified according to this line 
quantum number As an example, in Fig.|5]we show that the 
uniform or "ferroelectric" ll42ll dimer covering corresponds to 
the vacuum of lines. Number conservation implies that this 
state cannot be connected with any other ground state by a 
local process. The other "ferroelectric" dimer state (of oppo- 
site polarization) which is obtained by spatial inversion corre- 
sponds in the line representation to the state which is full of 
lines, and by the same argument is not connected by a local 
process to any other states in the ground manifold. Spatial 
inversion corresponds in the line representation to the line- 
antiline transformation (line conjugation). 

Not unexpectedly, the exact solution of the six-vertex 
model j42S] reveals that the ground states of greatest impor- 
tance (by which is meant those becoming increasingly sharp in 
the thermodynamic limit) are those containing, on a checker- 
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FIG. 5: (Color online.) Uniform dimer covering or "ferroelectric" 
state l43l . coffesponding to vacuum of lines. The opposing "ferro- 
electric" polarization obtained by spatial inversion corresponds to the 
state which is full of lines. 



board lattice of size L x L , precisely L/2 lines. This sector 
has the highest statistical (entropic) weight by virtue of the 
many ground states it possesses. Analogously, we will see in 
Sec. V that when the system is perturbed away from the Klein 
point at T = the state stabilized by the perturbation is gen- 
erally one with L/2 lines ll47ll . 



C. Critical correlations 

To expand upon the previous statements, we stress that all 
of the dimer coverings in the ground-state manifold are in- 
dividually eigenstates of the Hamiltonian. Thus there are 
no quantum fluctuations at the Klein point, and a zero- 
temperature transition through this point, for example as a 
function of K/J, is a discontinuous, first-order quantum 
phase transition j:7':]. The Klein point is a distinct type of crit- 
ical point which is essentially classical in nature with critical 
thermal fluctuations at all temperatures T > Q. 

It has long been known ll42ll that the square-lattice ice 
model formed by the centers of the crossed plaquettes in the 
checkerboard lattice exhibits critical power-law correlations, 
which have been investigated extensively in Refs. m. Ele- 
mentary entropic arguments combined with the ice condition 
lead to a simple dipolar analogy suggesting that the dimer cor- 
relations in d dimensions exhibit an decay, where r is 
the separation between dimers UtIi . Specifically, denoting 
by P = (±1,0,0) the first two six-vertex configurations of 
Fig. [3l and similarly by P = (0, ±1,0) and P = (0, 0, ±1) 
the other two pairs, the line representation applied for the py- 
rochlore lattice yields ITtIi 



(15) 



for large separation |r|, with i,j = 1, 2, 3 denoting the spa- 
tial components of P and r the unit vector in the direction of 
r; A is a constant. The correlation function for the checker- 
board lattice has a similar |r|^^ decay. When interpreted as 
flux lines, the six-vertex arrows (Fig.|3]l adhere at every ver- 
tex to a strict condition of no divergence. This correspon- 
dence, along with the entropic arguments, underlies the dipo- 



lar correlations. The zero-flux condition at every vertex leads 
to the long-ranged correlations (in fact with infinite correla- 
tion length) of the dimers. 

The ground-state manifold at K = Kc coiTesponds to 
the high-temperature limit of the ice model, in which every 
six-vertex configuration is realized with equal amplitude. As 
such ll42ll . the system is in the "disordered" critical phase de- 
scribed earlier where, in particular, the expectation value of 
any local six-vertex order parameter vanishes ll42ll . On the 
checkerboard (pyrochlore) lattice the different sets of states 
in the manifold remain characterized by the net number of 
vertical lines on each horizontal row (plane) of lattice sites 
(or conversely): this number is conserved and corresponds to 
a topologically invariant "flux" of lines. Because lines of ar- 
rows cannot terminate in the six-vertex representation the flux 
cannot change from one plane to the next in the pyrochlore 
ll42ll . We remark for clarity that for the checkerboard lattice 
the "horizontal" and "vertical" directions specified above for 
the dual lattice are rotated by 45*^ relative to the horizontal and 
vertical directions of the original lattice. Although the system 
does not possess local order it displays this elementary topo- 
logical order characterized by the topologically invariant flux 
ll48ll . The model discussed here forms a high-dimensional 
realization of Wen's string nets li48ll . On the checkerboard lat- 
tice, conservation of the line quantum number means that if 
ny is the number of vertical lines intercepting a horizontal 
span of fixed length at position Y on the y axis, the non-local 
correlation function for two strings at y = Fi and y = Y2, 
which is given by 



(16) 



with (f) arbitrary, is maximal, G(Yi, Y2) = 1. Because in gen- 
eral \G\ < 1, for unrestricted values of ny^ and nyj, here 
one has simply rty^ = ny^ . Similar results apply for the py- 
rochlore lattice, where ny^ denotes the number of lines inter- 
secting the plane Y ~ Yi. Equation ( fT6b presents a high- 
dimensional analog of the well-known fixed value of string 
correlation functions in 5 = 1 spin chains at the AKLT point 
ll49ll50ll . At large separations these uniform non-local corre- 
lations dwarf the local dimer-dimer correlations, which retain 
an algebraic decay (as r^"^). As a direct consequence (below), 
the system exhibits in general a long-ranged dimer order away 
from the Klein critical point. 



D. Low-energy excitations at the Klein point: spinons 

The energy cost of exciting one of the singlets of a given 
dimer configuration is proportional to J; as depicted in Fig.|6] 
this energy is not altered by changes in the separation of the 
two 5=1/2 objects forming the triplet state, which may thus 
be regarded as two deconfined spinon excitations. Single- 
spinon propagation can be described as a constrained random 
walk because the allowed paths may not increase the number 
of tetrahedra without a singlet dimer, but in contrast to the 
situation in Ref. |@| the paths are d-dimensional rather than 
being constrained to a one-dimensional path. 
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FIG. 6: Example of spinon propagation. The potential energy does 
not depend on the length of the string connecting the pair of spinons, 
which consequently are in a deconfined phase. In the line represen- 
tation the position of each spinon is indicated with a small dotted 
circle. The shaded square corresponds to a plaquette without singlets 
that we denote as a "defect". 



In this context it is instructive to consider the dimensional- 
ity of the allowed spinon paths for the checkerboard and py- 
rochlore Klein models, to examine its relation to the dimen- 
sionality dg of the symmetry group emerging in the ground 
state sector |@] . Here dg is the dimensionality of the minimal 
non-empty set of spins influenced by these symmetry opera- 
tions. Assuming that the ground-state sector does not trans- 
form trivially under these symmetry operations, the entropy 
Sg scales according to Sg ^ N'^"^'^, where dg ^ d — dg. 
The Klein model on the pyrochlore lattice exhibits a zero- 
dimensional {dg = 0) symmetry, meaning that there exist 
local operations which link different ground states (e.g. the 
operation depicted in Fig.|4]i. 

The presence of these symmetry operators enables spinon 
excitations to propagate with no energy cost in dg dimensions. 
Thus spinons on the pyrochlore and checkerboard lattices 
propagate freely in regions whose size scales with the vol- 
ume of the entire lattice, dg = d. For the square-lattice model 
of Ref. 120, d ~ 2 while dg = 1, as a consequence of which 
the spinons can propagate only in dg ^ d — dg ^ 1 dimen- 
sional regions (lines). This analysis provides a more explicit 
definition of the terms "complete" (pyrochlore and checker- 
board) and "partial dimensional reduction" (square lattice), of 
their origin in terms of system entropy (scaling respectively 
with the volume and with the perimeter), and of their conse- 
quences for spinon dynamics. 

As noted earlier, perturbation-theoretic calculations, in- 
cluding those illustrating the effects of spinons, are somewhat 
involved in the dimer basis due to the non-orthogonality of 
the dimer states. In the remainder of this section we comment 
only briefly and qualitatively on possible spinon exchange 
paths connecting different dimer states of the low-energy sec- 
tor; a detailed discussion of one systematic approach to quan- 
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FIG. 7: (Color online.) Spinon propagation in the "ferroelectric" 
or uniform state depicted in Fig. |5] Disruption of the dimer back- 
ground dictates that the line or string created between the propagat- 
ing spinons can move only in the direction indicated by the black 
arrows. This implies that the spinons are restricted to move within 
a cone defined by the horizontal and vertical directions. In the line 
representation, the position of each spinon is indicated with a small 
dotted circle. 



titative calculations in the dimer basis is deferred to Sec. V. 

While we have used the words "free propagation" of 
spinons to indicate that their motion occurs in all spatial di- 
mensions, spinon dynamics on the pyrochlore and checker- 
board lattices at the Klein point remain constrained by the 
effective ice rules obeyed by the dimer configurations in the 
ground-state manifold. These rules, upon which we remark 
further below, become most evident by considering the line 
representation of the six-vertex model (Fig. O. In the pres- 
ence of a single dimer excitation, the allowed lines specify the 
possible spinon motion. Against the background of the uni- 
form "ferroelectric" state (corresponding to the line vacuum. 
Fig. IS, where all horizontal and vertical arrows point in the 
same directions, the spinon lines which may be inserted are 
all chiral. The spinons are restricted to propagate only within 
a cone defined by two spatial quadrants (Fig.|7]i, and no closed 
spinon paths are possible in such states. The same is true for 
the opposing "ferroelectric" state, the state full of lines, in 
which the string created between the two spinons appears as 
an antiline. 

As suggested by the simple example of the previous para- 
graph, the spinon motion can be determined from the line 
rules. The trace left by two spinons which are created on a 
common tetrahedron and move away from each other corre- 
sponds, in the line representation, to a string which changes 
the state of the underlying link: the string creates a line seg- 
ment if the link was empty and annihilates the corresponding 
segment if it is already present. Following the rules of Fig.|3] 
if the line segment created by the string crosses a preexisting 
line, the resulting state on the corresponding tetrahedron is 
that indicated as 2 in Fig. [3] Once again we show for clarity in 
Figs.ISc) and (d) the line representations of the dimer/spinon 
states depicted in Figs. |6|a) and (b) for the checkerboard lat- 
tice. In the general case of non-uniform states with a finite 
density of lines, the spinon strings are no longer constrained 
to be chiral (Fig.|6]l. When the string of a propagating spinon 
collides with a six-vertex line representing the dimer config- 
uration, the lines (and hence the spinon path) are effectively 
reflected. 
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IV. THERMALLY DRIVEN DECONFINEMENT AWAY 
FROM THE KLEIN POINT 

As noted above, the dipolar correlations characterizing the 
deconfined Klein critical point are driven by classical fluctu- 
ations, whence the most appropriate designation is as a clas- 
sical r = 0+ critical point. Away from T = 0, the critical 
thermal fluctuations of the dimer coverings produce an en- 
tropic spinon-spinon interaction. This interaction may be de- 
termined by replacing each spinon with a static monomer in- 
serted in the dimer covering, reducing the problem to the clas- 
sical statistical mechanics of close-packed dimer coverings 
lf5lll52ll . This reduction makes use ofjhe fact that the singlet 
coverings are linearly independent |53f|. The critical dipolar 
fluctuations of the underlying dimer field produce an effec- 
tive "Coulomb" interaction between the spinons ifl?!], which 
in the pyrochlore lattices under consideration arises from the 
local conservation law of zero divergence. 

We will demonstrate the existence of the Coulomb phase 
for the checkerboard lattice; a similar analysis may be ap- 
plied for the 3d pyrochlore case. The derivation of the ef- 
fective Coulomb interaction begins by assigning a charge to 
each spinon, for which it is convenient to expand the six- 
vertex representation in order to include states containing two 
spinons. It is important to note that spinons are created in pairs 
(singlet-triplet excitation), and each pair generates a single 
crossed plaquette of higher energy, because is possesses no 
singlets, which will be called a "defect". When the spinons 
propagate away from each other they are also separated from 
the defect (Fig.|6]l, and preservation of the ice-rules requires 
that the number of defects does not increase. We will show 
below that the effective charge of the defect is opposite 
in sign to the effective charge Qs of the spinons, with magni- 
tude Qd = — 2(7s. [Superficially, this overall neutrality of three 
quasiparticles forms a geometric analog of the quark content 
of the neutron, whose basic quark structure is \udd) with the 
charges of the u and d quarks being respectively (2e/3) and 
(-e/3).] 

To represent states with two spinons, we will include the 
four additional vertex configurations depicted in Fig. [8] Note 
that these configurations have non-zero divergence, in con- 
trast to the six-vertex states (Fig. O used to represent the 
dimer coverings (which are the vacuum of charge). According 
to the previous rules, the defect is represented by four outgo- 
ing arrows [Fig. [Sja)]. For the checkerboard geometry it is 
convenient to introduce the two Neel sublattices A and B of 
the underlying square lattice; a number of equivalent defini- 
tions may be found for a two-sublattice decomposition of the 
3d pyrochlore. A one-to-one mapping between states con- 
taining spinons and the vertex configurations may be found 
by restricting the spinons to one of the sublattices (for exam- 
ple A). A complete classification of all spinon configurations 
may be obtained by introducing an additional Ising-like quan- 
tum number which specifies on which of the two sublattices 
our plaquettes are located. Such an additional tabulation is 
not necessary for the examination of the long distance inter- 
action. We may rely instead only on the consequences of the 
vertex representation illustrated in Fig.[8]for plaquettes of one 
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FIG. 8: (Color online.) (a) Vertex representation of a crossed plaque- 
tte without any singlet (defect), (b) Vertex representation for spinons 
that are located in one sublattice of the checkerboard lattice. 



sublattice, which are sufficient for determining the asymptotic 
(long-distance) form of the interspinon interaction. The fact 
that spinons may propagate from one sublattice to the other 
due to the existence of diagonal bonds in dimer coverings of 
the checkerboard lattice is nevertheless crucial. This is not the 
case for quantum dimer models on bipartite lattices, whose 
dimers always link the two different sublattices. We remind 
the reader that on a bipartite lattice each spinon may be as- 
signed a unique effective charge given by the sublattice parity 
of the spinon location (usually H-1 on one sublattice, and -1 on 
the other). This simple assignment is invalid on non-bipartite 
lattices such as the pyrochlore, where the additional range of 
possible spinon configurations necessitates an essential exten- 
sion of the definition of the effective spinon charge. 

Because the spinons occupy sites of the original lattice, 
their position will be denoted by one arrow on the dual lat- 
tice. Specifically, we take this as the ingoing arrow coming 
from the lower left corner of a plaquette with only one outgo- 
ing arrow [Fig.[8tb)]. If one considers vertex states such that 
for each vertex with maximal divergence [Fig.[8ta)] there are 
two with only one outgoing arrow [Fig.Ob)], the mapping be- 
tween the two representations is one-to-one (under the con- 
dition that all the spinons occupy the same sublattice). The 
Coulomb interaction and the effective charge of the spinons 
emerge from the standard procedure used to map the ice prob- 
lem onto a compact U(l) gauge theory llsill . Here we follow 
the derivation and notation of Ref. 15511 to describe this proce- 
dure and derive the effective charges of defects and spinons. 

The local zero-divergence condition defining the dimer 
coverings on the checkerboard and pyrochlore lattices implies 
that these states can be considered as the physical Hilbert 
space of a compact U{1) gauge theory. The local U{1) gauge 
invariance becomes more explicit if one defines pseudospin- 
1/2 operators Sf on each site i of the lattice. The dual lat- 
tice formed by the centers of the tetrahedrons (vertices) can 
be divided into two sublattices A and B. The site coordinate 
i of the original lattice denotes the links of the dual lattice. 
5f = 1 /2 if the arrow on link i points from the A to the B 
sublattice, iSf = —1/2 if the arrow points from B to A. The 
zero-divergence condition can be expressed in terms of these 
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variables as 
and the operators 

6^-^° ^-eo 5r (18) 

are the local U(l) symmetry transformations associated with 
this conservation law. We note that the conservation law is 
a property of the Hilbert space and is therefore present for 
any theory defined on this space. An equivalent constrained 
low-energy Hilbert space was derived for an XXZ model on 
a pyrochlore lattice in Ref. lEoll . 

We present in some detail the consequences of the zero- 
divergence condition for the checkerboard lattice (d = 2), for 
which the notation is slightly more compact; the results we 
obtain are unchanged for the pyrochlore lattice {d = 3). The 
U{1) gauge structure of the Hilbert space indicates that one 
may gain further insight by identifying the variables playing 
the role of the electric field and of the vector potential in elec- 
trodynamics. For this purpose, it is convenient to introduce 
quantum rotor variables and the coordinate r = nei + me2 to 
span the dual lattice (the lattice of centers of the tetrahedrons). 
The primitive vectors ei and 62 define the two perpendicular 
directions on the dual lattice. On each link of this lattice we 
define an angular-momentum variable /^(r) = ±1/2, where 
/i(r) indicates the oriented link from r to a neighboring dual- 
lattice site r + e^. The state ^^(r) = 1/2 ([^(r) = -1/2) 
corresponds to an arrow parallel (anti-parallel) to the relative 
vector e^j. With this definition, /^(r) changes sign under in- 
version: /_p(r) = — ^/j(r) (— /i(r) indicates the oriented link 
from r to r — e^). For convenience we introduce the variable 
l^{r) = l^{r) + 1/2, which takes the integer values or 1. 
We employ the discrete gradient 

V^/,(r) =;,(r)-Z,(r-e^) (19) 

to reexpress the ice rules as an explicit zero-divergence con- 
dition on the vector field l^{r), 

V^/^(r) = 0. (20) 

Equation (|20] i corresponds to Gauss' law for electrodynam- 
ics in the absence of external charges, V • E = 0. Ver- 
tex configurations violating the ice rules, such as those de- 
picted in Fig. [8] should then carry an "effective charge" 

Q = E^=i,2^a'^m(i")- Here Q = E,,=i,2^p'p(i") = 
2 for the configurations of Fig. |8ja) (defects) and Q = 

J2p.=i 2'^ fj-^fj-i^) = ^1 ^'^^ three configurations of 
Fig. |8jb) (spinons): the spinons carry a positive unit charge 
while the defects carry two units of negative charge. The 
"electric field" is then defined as the operator E^(r) with 
eigenstates |^;^(r)) and eigenvalues Z,i(r). The canonically 
conjugate operator 6,y(r) satisfying 

[e,(r'),E^(r)] =iJr,r'V>^ (21) 



plays the role of the vector potential. Because in this section 
we focus only on the high-temperature (T 00) fixed point 
whose properties are determined solely by the gauge struc- 
ture of the Hilbert space, we will not discuss any particular 
quantum gauge theory. The above derivation is used only to 
identify which variables play the role of the electric field and 
the charge. 

The Klein point at T = 0+ corresponds to the infinite- 
temperature fixed point. Such a point is characterized by a 
deconfined two-component (spinons and defects) Coulomb 
gas. The spinons exist only as gapped excitations (of gap 
A oc J), and the spinon-spinon interaction is then logarith- 
mic, V{r) = jJ2j<k Ijik In i^jk/a), for d = 2 and a power 
law, V{r) = 7X]j<fe QjQk/rjk, for d ~ 3. The spinon-defect 
interaction has the opposite sign (attractive) and its absolute 
value is two times greater Although a spinon-defect pair is 
critically confined in the d = 2 case by the logarithmic inter- 
action, the prefactor 7 is generally sufficiently small that there 
is a divergence in the expectation value of {r"^), which im- 
plies a divergent "dielectric constant" lf56ll . At finite tempera- 
tures there is a finite concentration p ^ e~^^ [(3 = (fceT)"^] 
of thermally excited spinons in equilibrium, and the effective 
spinon-spinon or spinon-defect interaction is screened, de- 
caying exponentially for distances r ^ However, the 
critical behavior remains observable for r <C l/p^/''. 

For general K ^ Kc and \K — Kc\ ^ J, the energy differ- 
ence between the singlet and triplet states on each tetrahedron 
is much smaller than the gap to the quintet state (S ~ 2). 
IstIi It is clear that for e = \K — Kc\ ^ J the degeneracy 
between the singlet and triplet states is effectively restored at 
temperatures T such that 

e < fceT < J. (22) 

In this case, only the manifold of singlet dimer coverings 
found as the ground states at K = Kc is accessible. Thus 
the effect of a finite temperature on a system in the proximity 
of a Klein point is the effective restoration of extensive con- 
figurational entropy, critical dipolar correlations, topological 
order, and spinon deconfinement. 

We comment that the appearance of deconfinement at high 
temperatures, also in systems far more general than those con- 
sidered here, is hardly surprising. We repeat well known 
facts: various ordered phases generally melt via the appear- 
ance of topological defects; these defects become deconfined 
at the melting temperature; the origin of deconfinement is an 
energy-entropy balance; entropic effects modify the effective 
interactions between the defects, leading to a vanishing inter- 
action at the onset of deconfinement. A careful treatment of 
entropic effects often involves detailed contour-counting ar- 
guments ifssll . The effective force between topological defects 
in a wide range of systems, including vortices in supercon- 
ductors and dislocations and disclinations in elastic media, is 
universally expected to be algebraic, as may be seen directly 
from duality arguments lf59l l60l I61I1 . In the current context 
the topological defects correspond to monomers or spinons in 
the dimer coverings and the Coulomb field lines to the many 
possible trajectories for the motion of a defect. Pure Coulomb 
interactions vanish at large separations in dimensions d > 2, 
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suggesting the possibility of high-dimensional fractionaliza- 
tion. The peculiar aspect of the pyrochlore systems under 
consideration is that the critical point separating the low- 
temperature confined state from the high-temperature decon- 
fined phase may be driven to T = 0+ by optimizing the frus- 
tration in H. This is the singular characteristic of the Klein 
point: the exponentially large ground-state sector is defined 
solely by the local zero-divergence constraint. The system is 
then driven only by entropic fluctuations, which in combina- 
tion with the local constraint give rise to the classical critical 
behavior 



V. ZERO-TEMPERATURE STATES AWAY FROM THE 
KLEIN POINT 




FIG. 9: (Color online.) The "antiferroelectric" state. Dots and 
crosses at the centers of the plaquettes of the dual lattice indicate 
the "flux" associated with the arrow circulation on the correspond- 
ing plaquette. The two antiferroelectric states are the only ones for 
which each plaquette of the dual lattice has a net chirality. 



A. Local perturbations in tlie non-ortiiogonal dimer basis 

In this section we examine explicitly the possible ground 
states away from the Klein point. To this end, we begin by 
developing a systematic procedure for evaluating "matrix el- 
ements", {4'a\0\tpb), within the non-orthogonal basis of sin- 
glet dimer coverings which constitute the low-energy sector 
in the vicinity of the Klein point. We note that because of 
the non-orthogonality a standard (bra-operator-ket) expecta- 
tion value computed between two states in the manifold is 
not a true matrix element, and is referred to henceforth as 
a "bracket". We will illustrate the calculational procedure 
and results for the checkerboard lattice, and state only that 
their generalization to the pyrochlore is straightforward (if 
lengthy). 

The evaluation procedure consists of elementary rules 
which are applied to the graph associated with the bracket to 
be evaluated, and it is a generalization of the standard method 
used to compute the overlap between two singlet dimer cov- 
erings on bipartite lattices \62:, ,63, 64, ^]. We first extend 
this method to compute the overlap between non-bipartite sin- 
glet coverings, because on the pyrochlore system the dimers 
are not necessarily formed from sites on different sublattices. 
Secondly and more importantly, we provide straightforward 
rules for computing brackets of spin-product operators of 
Heisenberg type. Si ■ Sj. Brackets of products of Heisenberg 
spin operators, as in the second term of Eq. ([T]i, can be evalu- 
ated by a simple extension of these rules. 

In this section, we will consider a perturbation which is pro- 
portional to the first term of Eq. a Heisenberg interaction 
between any pair of sites in the same crossed plaquette 

iJi = AJ ^ 5f • Sf, (23) 

such that after the perturbation K ^ 4J/5. For the Klein 
model on the square lattice it is transparent that all states 
within the ground-state basis remain degenerate to first order 
in the perturbation Hi. This statement follows from two ob- 
servations: the diagonal brackets, (V'a|-ffi|'0a)> the same 
for any dimer state \ipa)', the off-diagonal terms, (V'a|C|'0b)' 
with a ^ b, are zero for any local operator O, in particular for 



O = Hi. Thus the degeneracy is lifted only to second order 
in Hi on the square lattice. The situation is completely differ- 
ent for the same model on the pyrochlore and checkerboard 
lattices: because in these cases the ground-state degeneracy 
is extensive at the Klein point, two different ground states can 
be connected by local operators. As a consequence the off- 
diagonal brackets do not in general vanish, and the degeneracy 
is lifted to first order in the perturbation. 

The usual approach to this problem is to neglect the non- 
orthogonality of the singlet dimer coverings and propose a 
minimal quantum dimer model which includes only processes 
of the Rokhsar-Kivelson (RK) type (6^ . We note, however, 
that this procedure does not lead to the RK model in the py- 
rochlore geometries under consideration because, although 
the Hamiltonian is the same (a U(l) gauge magnet ll66ll ). the 
zero-divergence or gauge condition which defines the "phys- 
ical" Hilbert spaces is different. In what follows we depart 
from this conventional treatment, and show that it is indeed 
incorrect for the case of interest. The non-orthogonality of 
the dimer basis introduces subtleties not directly evident from 
the conventional approach, as a consequence of which we will 
show that the valence-bond crystal (VBC) orderings chosen 
by tunneling events differ from the ones expected for a theory 
which includes only processes of the RK type. 

In terms of the arrows (the six-vertex representation) 
shown in Fig. |9l the antiferroelectric state becomes a stag- 
gered flux phase in which each square plaquette of the dual 
lattice has a well-defined chirality (all of its arrows circulate 
in the same way) and the direction of circulation alternates 
from clockwise to counterclockwise. With this starting state 
one may invert the sign of the current (direction of the arrows) 
around any plaquette. This is the elementary process, or cor- 
ner flip in the line representation, depicted in Fig. |4|b). The 
antiferroelectric state allows for the largest number of these 
elementary processes, and each can be interpreted as the cre- 
ation of a localized defect on an uncrossed plaquette of the 
antiferroelectric background. The tunneling between states 
with defects in different positions lowers the energy of the 
system. While the creation of one of these defects is indeed 
the elementary RK process normally invoked as the most rel- 
evant one in these cases ll39ll . we will show that its amplitude 
is exactly equal to zero for the perturbation Hi. 



10 



February 6, 2008 



submitted to Physical Review B 



To first order in the perturbation the state of lowest energy is 
the linear combination of singlet dimer coverings which min- 
imizes Hi . Up to an irrelevant constant, the Hamiltonian Hi 
can be rewritten in the more convenient form 



Hi = lAj ^ p.r„ 



(24) 



where P," is the permutation operator which interchanges the 
spins on sites i and j of plaquette a. The diagonal brack- 
ets are manifestly the same for all singlet dimer coverings. 
Specifically, on the checkerboard lattice and for any system 
with PBCs, 



Ed = {'lPa\Hi\ija) = iNAJ, 



(25) 



with N the number of lattice sites and \tjja) denoting any 
dimer covering satisfying the Klein constraint (all crossed 
plaquettes contain one dimer). As noted in Sec. Ill, sys- 
tems with OBCs possess additional ground states containing 
crossed plaquettes with two dimers or (one or two) spinous, 
and these would have a different diagonal energy but are not 
thermodynamically relevant. We focus henceforth on ground 
states satisfying PBCs. The degeneracy of dimer coverings 
is lifted only by superposing different allowed coverings with 
finite amplitudes, i.e. by considering states of the form 



with the normalization condition 



5Ilcfcl' + I]c;cfe(V'd^fe) = 1. 

k l^k 



(26) 



(27) 



B. Loop Structure and Rules 

We proceed by listing a complete set of rules for the evalua- 
tion of brackets in the dimer basis. These general, purely topo- 
logical rules extend early and pioneering ideas concerning 
resonating valence-bond states, and initial results obtained 
by using these concepts for square lattices ||62[ l63l l64i l65ll . 
The normalization of and computation of {ip\Hi\ijj) re- 
quire evaluating the overlaps (?Aa | V'fo) and brackets of the form 
{'4'a\Pij\'4'b), with \ipa) and |-0b) any states in the dimer basis. 
To compute the overlap we first assign a given orientation to 
each singlet dimer, indicated with arrows pointing upwards 
for vertical dimers and to the right for diagonal and horizontal 
bonds (Fig. [TOl i. [We remark that "vertical" and "horizontal" 
refer here to the real lattice, and not to the rotated lattice used 
in the discussion of the six-vertex and line representations in 
Sec. 111.] With this convention, a superposition of dimer con- 
figurations corresponding to \ipa) and \ipb), with the differing 
bonds marked as in Fig. fTOl c) (bonds which do not differ give 
trivial loops of length two, introducing factors of unity), forms 
a graph G"''. In general, this graph may consist of more than 
one disconnected loop F"''. 
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FIG. 10: (Color online.) Example (c) of the graph G"'' which results 
from superposing two different dimer coverings [(a) and (b)]. These 
dimer coverings correspond to an antifen'oelectric background with a 
single elementary process operating on different (nearest-neighbor) 
uncrossed plaquettes. 



Let L{Tj'') be the length of the loop F"^ (the number of 
singlet dimers forming the loop), iVr(G°'') be the number of 
disconnected loops in graph G°^, and ric be the number of 
arrows in the graph which circulate clockwise in their loop. 
With these definitions the overlap {ipa\tpb) may be expressed 
as 



where 



L(G"^) 



(28) 



(29) 



j = l,Afr(G»'>) 



is the total length of the graph. The determination of the non- 
trivial brackets {^a\Pij\^b) proceeds by examining how the 
length and dimer orientation of the graph G"^ are affected 
when the permutation operator Pij is applied to the state \ipb)- 
The resulting graph will be denoted as Gfj. The application 
of Pij may alter the loop geometry in only a restricted set of 
simple ways, of which we next provide a systematic list. 

[i] If i,j e G"** and they belong to the same loop F, there 
are two possibilities: 

(a) when the distance (number of dimers) on the loop 
between sites i and j is an odd number, the effect 
of Pij is to cut the loop into two segments and 
to reconnect them while exchanging the ends of 
one of the segments [Figs. fTTTa) and (b)]. The 
length of the loop under consideration is unaltered 
in this process. The only effective modification is 
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the reversal of an odd number of arrows, causing 
a sign-change relative to (V'alV'fc)' 



{'ipa\Pij\lpb) = -gab- 



(30) 



(b) when the distance between sites i and j is even, 
Pij divides the loop F into two smaller loops, 
Fi and F2, preserving the total length, £(F) ~ 
L(Fi) + L{T2) [Figs. [TTTc) and (d)]. The circu- 
lation of the arrows is unaltered, and thus from 
Eq.(l28Tl the net change relative to graph G"^ is a 
supplementary factor of 2 due to the increase in 
the number of loops, Nr{Gf^) = iVr(G"'') + 1, 
whence 



{'ipa\Pij\'ipb) = '^gab- 



(31) 



[ii] If only one of the two sites (i or j) belongs to graph G"'', 
the effect of Pij is to increase the length of the graph by 
two Unks, /.(G^/) i(G"^) + 2. There is no change 
of sign because the arrows of the two additional bonds 
circulate in opposite directions relative to their point of 
contact. From Eq. (|28] | the net result relative to graph 
G"** is a factor of 1/2 due to the length increase, and so 



\gab- 



(32) 



[iii] If neither i nor j belong to graph G"'', there are two 
possible cases: 

(a) if there is a singlet bond connecting the sites i and 



(V'alPylV'b) = -gab] 



(33) 



(b) if the sites i and j are not connected by a singlet 
bond. 



mPjH'b) 



\gab- 



(34) 



[iv] Finally, if i, j G G°^ but the sites belong to differ- 
ent loops Vi and Fj, the effect of is to fuse both 
loops into a single loop F^ . Because the total length, 
L(Fij) = L{Vi) + L{Vj), and the sign are preserved 
under this operation, the net result is a factor 1/2 which 
arises from the decrease by one in the number of loops. 



Nt{G1°) = NriC") - 1, leading to 



{MPMb) 



\gab- 



(35) 



These rules allow a very straightforward derivation of 
Eq. ( |25T l: because G"" is the null graph, the effect of cor- 
responds to case inil. For any dimer covering there are N /2 
bonds of type liiial and hN/2 bonds of type liiibl . and adding 
the two contributions yields 

Ed - mHl\i^a) - (-gaa^ + • (36) 
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FIG. 11: (Color online.) Effect of the permutation operator Pij on 
graph G"'' (Fig. 1101 ) when sites i and j belong to the same loop of 



G" and their separation is odd [case (a) 
(d)]. 



(b)] or even [case (c) 



From gaa = 1 [see Eq. ( |28] )1 we recover Eq. ( |25T l. 

With the same rules one may compute the energy of the 
most general state within the Klein-point ground-state mani- 
fold whose expression is given by Eq. ( |26l ). 



l^k 



/3i(G"^) 
gik — ^ 



IV 



lk\ 



Wlk 



Here Wit 
pairs 

G"=, while V 



(37) 

n'l'^, where n'J; (n'^) is the number of 
which are of type ITa) (lEbl) relative to the graph 



Ik 



^ is the total number of pairs 
of type [i]. The reason why only cases of type [i] need be 
considered in Eq. dJTl ) is that terms of types Inl, HiibL and 
iHvl all give the same contribution, gab/"^- 



C. Energy Calculations 

Next we apply the rules of Subsec. VB to the checkerboard 
lattice. Equation ( l37b allows one to compute the energy of 
any variational state. Its second term corresponds to the en- 
ergy change 5E due to the linear combination of different sin- 
glet dimer coverings. The magnitude of this change therefore 
reflects the degree of "mixing" between different dimer cov- 
erings induced by Hi . To develop an intuitive understanding 
we begin by considering the linear combination of two dimer 
coverings, a natural choice being to superpose two coverings 
whose dimer configurations differ minimally (such as the two 
states depicted in Fig. [Hi, i.e. by the RK or corner-flip process 
discussed above. We remind the reader that the six-vertex 
representation of this process entails a reversal of the arrows 
on a single plaquette of the dual lattice with a well-defined 
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FIG. 12: (Color online.) Graph G°'^ resulting from the superposition 
of two dimer coverings which differ minimally. The physical process 
connecting these is the RK-type or corner-flip process represented 
in Fig. m 
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FIG. 13: (Color online.) Graph G° resulting from the superposition 
of two states connected by inverting the arrows of two next-neighbor 
plaquettes of the dual lattice which have opposing, well-defined chi- 
ralities. The two plaquettes are highlighted in panels (a) and (b). 



chirality [Figs.[T2ja) and (b)]. If |-0o) and \'ipd) are any such 
pair of dimer coverings which differ minimally, their overlap 
is 

god = (V'olV'd) = 2x2-4=1, (38) 

because L(G°'') = 8 and = 4 [Fig. [He)]. To obtain the 
energy of the normalized state 

\i^od)=a\tlJo)+PWd), (39) 

one replaces the coefficients of Eq. dJTj i with the values dic- 
tated by the above rules. The number of pairs (z, j) with both 
sites on the graph is |V°''| = 12, composed of 8 pairs of type 
Hal and 4 of type llbl, as a result of which W°'^ = 0. Replac- 
ing these quantities in Eq. ( [37] i gives 

{^bod\Hi\i'od) = Ed. (40) 

The fact that the energy is unchanged {SE = 0) by taking 
the linear combination ( [39] l indicates that the elementary RK- 
type process is not generated by Hi . We note that all pairs of 
sites {i,j) in the only loop of the graph G"** [Fig. [T2l c)l are 
separated by a distance on the loop no greater than 2; loops 
satisfying this condition will be defined as simple. As an ex- 
ample, the loop depicted in Fig. [TOl c) is not simple because 
it contains site pairs separated by five and six bonds. 

In general, 6E vanishes for any pair of dimer configurations, 
iV'a) and \ipb)^ such that all loops in the graph G"'' are simple. 
This property follows immediately from Eq. ( |37] |: if all loops 
are simple one has that n"^ = 2nfl^ = L(G°^), implying 
Wab = and |V°^| = 3L(G'''')/2. 

Because Hi does not connect two minimally different 
dimer configurations [L(G°'') = 8], it is necessary to seek 



the minimal relevant process, meaning that with the minimum 
value of L(G°^), which is generated by Hi. The next possible 
value is L(G''^) = 12 (no graphs exist with L(G"'') = 10), an 
example of such a physical process being that connecting the 
states (a) and (b) in Fig. [10] In the six-vertex representation 
this process is permitted when all the arrows of a rectangle 
formed by two adjacent plaquettes of the dual lattice have the 
same chirality; this rectangle is highlighted in Figs. [TOf a) and 
(b). The process consists of inverting the orientation of the 
six arrows, or the flux of the rectangle. For any pair of dimer 
coverings ji/'o) and ji/'i) which differ by this process, 

ffoi = (V'olV'i) = -2 X 2-6 = (41) 

and for the normaUzed linear combination \^qi) = co|V'o) + 
ci the energy change is 

(V'oiliJil^oi) = + ^coCiAJ. (42) 

The minimum energy is obtained for |co| = |ci|, which to- 
gether with the normalization condition gives cq = ci = 
l/v/2(l+.goi)forAJ < Oandco = -ci = 1/^2(1 - 301) 
for A J > 0. The overlap is \goi\ = 1, giving 

I Co I = |ci| ^ l/\/2, as expected for the mixing of two or- 
thogonal states with the same diagonal energy. This remains 
true for any pair of states connected by Hi because jgoi | is an 
upper bound for the absolute value of their overlap. 

In the next process to consider, the arrow directions on the 
perimeters of two next-neighbor plaquettes of the dual lattice 
are inverted if both plaquettes have opposing, well-defined 
chiralities, as indicated in the examples of Figs. \T3\a) and 
(b). In this case we denote the corresponding dimer cover- 
ings by iV-'o) and {^2), and the length of the associated graph 
is L(G°^) = 14. The overlap is 502 = 1/64 and the en- 
ergy shift due to the perturbation for the linear combination 
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FIG. 14: (Color online.) Graph G"'' resulting from the superposition 
of two states connected by inverting the arrows of two next-neighbor 
plaquettes of the dual lattice which have the same well-defined chi- 
rality. The two plaquettes are highlighted in panels (a) and (b). 



|V''02) = ColV'o) + C2\i>2) is 



{ip02\Hi\ipo2) = Ed - Y|gCoC2AJ, 



(43) 



with Co = ±C2 = 1/^/2(1 ±.902) for A J = ±|AJ|. 

The other graph G"'' of length L{G°'^) = 14 corresponds 
to a process which is analogous to the previous case but with 
next-neighbor plaquettes of the dual lattice having the same 
well-defined chirality (Fig. \T4\ . The energy change SE as- 
sociated with this process is exactly that obtained for |V'02)- 
Any other process generated by Hi has an associated graph of 
length L{G'^^) > 14, and because the amplitude of these pro- 
cesses falls in proportion to ?i°^2~^('^°''^/^ (72"^^ is the number 
of pairs {i, j) on the same loop separated by a distance greater 
than 2), these may safely be cut off beyond a certain length 



D. Variational Wave Function 

In order to propose a wave function which takes ad- 
vantage of the three processes giving non-zero contributions 
for L{G°-^) < 14, we introduce the bosonic operators Ag 
and Ag. To define these operators we classify the uncrossed 
plaquettes of the checkerboard lattice according to the sublat- 
tice A OT B (Sec. IV) of the dual square lattice [Fig. [T5la)l. 
In the six-vertex representation of the antiferroelectric state 
(Fig. |9]), these two sublattices have opposite chirality, corre- 
sponding to a staggered flux phase. If the uncrossed plaquette 
a has a well-defined chirality, the operator A^ raises the flux 
from negative to positive if a E A and lowers it from positive 
to negative if a E B, i.e. Ag changes the direction of arrow 
circulation from counterclockwise to clockwise if a G A and 



conversely if a E B. For states of no well-defined chirality 
on plaquette a, and when the flux is positive (negative) for 
a E A (B), Ag becomes the null operator, A^ = 0. 

Minimizing the energy of a variational state is equivalent to 
maximizing the number of allowed tunneling processes with 
a negative sign. The appropriate dimer background for this 
optimization is provided by the antiferroelectric state j-^afc) 
of Fig. |9j only for this state can the elementary process of 
Fig. [12] be applied to every uncrossed plaquette. Moreover, 
the three processes which give finite contributions to SE for 
L{G°'^) < 14 can be applied to any pair of nearest-neighbor 
(Fig.fTOl) and next-neighbor (Figs.[T3]and[T4li plaquettes of the 
dual lattice when the antiferroelectric background is suitably 
populated with the local defects (bosons) created by A^. For 
positive A J, these processes lead to a uniform energy gain 
(SE < 0) when the bosons are created with opposite phases 
on the sublattices A and B [Fig. [T5l a)l. The considerations 
for A J < mirror those for positive A J with a change of the 
decomposition in two sublattices: for negative A J the uni- 
form energy gain occurs when the bosons are created with 
opposite phases on the sublattices C and D [Fig. [T5l b)l. Em- 
ploying the definitions and observations introduced earlier in 
this section, the simplest wave functions for approximating 
the ground state in the presence of Hi are 



IV-) = t;^ n (1 + ^a) n (1 - ^a)l^afe), (44) 



where |i/'±) correspond to AJ = ±|AJ|. The coefficients 
u and V are determined from the normalization of the corre- 
sponding wave functions [Eq. (|27]il. Fig. [T5] illustrates the 
sublattice decompositions of the checkerboard geometry for 
the two cases: {A,B) for AJ > (|V^+)) and (C, D) for 
A J < (|V^_)). A similar decomposition into plaquettes C", 
D' is obtained by a 7r/2 rotation of Fig. fTSl b). 

The wave functions and |V'-) represent linear com- 
binations including all possible two-plaquette resonance pro- 
cesses, the phase relations of which are specified by Fig. [15] 
Although these are superpositions of many states from the 
manifold of dimer coverings, the number of participating cov- 
erings remains a very small fraction of the total, and the se- 
lected linear combinations are VBC states. In the pseudospin 
T = 1/2 space generated by the two possible states of each 
uncrossed plaquette of \ipi^ic), |V'+) and |i/'_) are analogous 
to the classical, magnetically ordered Neel [q = (tt, tt)] and 
collinear [q = (0, tt) or g = (tt, O)] states polarized in the x- 
direction. In this space j^/'afc) corresponds to one of the fully 
polarized ferromagnetic states (such as r? = —1/2 for all a). 
However, and \ t/j-) are not exactly the classical Neel and 
collinear states because two nearest-neighbor uncrossed pla- 
quettes cannot be "flipped" independently of each other, i.e. 



afc/ 







(45) 



if a and a' are nearest-neighbor uncrossed plaquettes. In the 
pseudospin representation the two plaquettes cannot be simul- 
taneously in the = 1/2 state, as a result of which these 
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FIG. 15: (Color online.) (a) Neel and (b) coUinear alternation of 
uncrossed plaquettes on two sublattices. (a) Decomposition of un- 
crossed plaquettes of the checkerboard system into two sublattices A 
(crosses) and B (circles), which dictates the form of the variational 
wave function J44t for A J > 0; (b) decomposition into sublattices 
C (crosses) and D (circles) which yields the form of the low-energy 
wave function for A J < 0. 



states are projected out of the superpositions corresponding 
to the classical Neel and collinear states to obtain and 
I ■(/'-)• We stress that the symmetries of and | -!/;_) for 
A J = ±1 AJ| are quite different. 

The energies of the wave functions and \t/j^) obtained 
from Eq. ( |37] | are 



E+ ^ {4>+\H,\^P, 



= Ed-^iL'N/^J 



Ed 



'^NAJ. 



(46) 



The three processes represented in Figs.[T0l[T3j and[T4]give an 
extensive negative contribution to E^ . The lack of frustration 
among these different processes indicates that \^+) and | -!/;_) 
are representative of the VBC orderings which are stabilized 
respectively for AJ > and AJ < 0. We draw particu- 
lar attention to the fact that in these valence-bond orderings 
the four dimers around each uncrossed plaquette resonate be- 
tween the two configurations shown in Fig. [12] The relative 
phase is positive on one sublattice (A for A J > and C for 
A J < 0) and negative for the other {B for A J > and D for 
A J < 0). Note that the RK process would lead to a different 
ground state in which this relative phase is the same for all 
the plaquettes of the dual lattice. In the pseudospin basis this 
corresponds to the ferromagnetic state [q = (0, 0)] with polar- 
ization in the x-direction for A J < and the — x-direction 
for A J > 0. 

With regard to possible four-spin perturbation interactions, 
we remark that the dominant type of such a perturbation 
would be that contained in the second term of Eq. ([Hi. How- 
ever, the effect of altering K is trivially identical to the per- 
turbation of Eq. ( l23T l, and the results derived throughout this 
section may be applied directly to Eq. ([D with a sign inversion 
of AK ={K- Kc) relative to A J. 

We now summarize the extension of our results to the py- 
rochlore lattice. Here any given dimer covering leads to pre- 
cisely the same, uniform energy shift (i25T l when a supple- 
mentary nearest-neighbor Heisenberg interaction is inserted. 
When superpositions of dimer coverings are considered, the 
elementary graph (or analog of the RK process) circumscribes 



Dilute plasma 
of spinons 




Kc K 

FIG. 16: Low-temperature phase diagram around the Klein point 
for checkerboard and pyrochlore spin systems. 



a hexagon formed by six neighboring tetrahedra. This elemen- 
tary hexagonal loop replaces the sides of the uncrossed pla- 
quettes on the checkerboard lattice. As for the checkerboard, 
the graph corresponding to the elementary process is simple, 
implying that this process is not generated by the perturbation 
Hi . Once again it is necessary to consider larger loops to find 
the valence-bond ordering stabilized by the perturbation. 

As stated above, the local zero-divergence constraint defin- 
ing the low-energy subspace around the Klein point implies 
that the low-energy theory induced by a general perturbation 
is always a compact U(l) gauge theory. Following the argu- 
ment of Polyakov, the phase stabilized by a general perturba- 
tion is always confining at T = 0. For the spin models we 
consider, this confining phase consists of a VBC configura- 
tion. The perturbation Hi treated in this subsection is a par- 
ticular example of this general observation, and also perhaps 
the most physically relevant. The effect of a general perturba- 
tion in the space of interaction parameters is illustrated in the 
schematic phase diagram of Fig.[T6l 

We conclude this section with a summary: we have in- 
troduced a loop-based approach to account for the physical 
processes in short-ranged quantum dimer systems with only 
local interactions, which is powerful and general. We have 
used this procedure to compute energies and variational wave 
functions in the non-orthogonal basis of dimer coverings on 
the checkerboard lattice. These considerations include all of 
the physically relevant perturbations of the Klein Hamilto- 
nian, as discussed in Sec. II, either explicitly or by direct anal- 
ogy. We conclude that the general case is one in which VBC 
phases, which are gapped and confined, are stable away from 
the Klein point. 



VI. EXACT CRITICALITY IN A FINITE-TEMPERATURE 
REGION 

At issue for the physical properties of a system is rather 
less the behavior at T = than that at finite temperature. In 
the preceding section we devised a systematic, purely topo- 
logical, diagrammatic method to compute the energies of all 
variational wave functions in the dimer basis. This procedure 
is exact for the degenerate manifold of basis states, including 
most importantly when this is non-orthogonal, and as shown 
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in Sec. V presents a tractable problem at zero temperature. 
At finite temperatures the need to evaluate brackets between 
all possible states in the Hilbert space presents a more seri- 
ous problem. In this section we discuss the extent to which 
progress is possible in this regime, by first illustrating how to 
establish the existence of a finite-temperature region of crit- 
ical behavior for general pertubations, and then consider a 
simple, anisotropic perturbation to (a) demonstrate explicitly 
the existence of this critical region and (b) derive its approx- 
imate phase boundaries. In the remainder of this section we 
set /cb = 1 and for rigor assume implicitly that J/T oo. 
J is the gap scale separating the ground state energies from 
the lowest excited states. T he p resence of a spectral gap in a 
Klein model was shown in 0911 . 



A. Finite-temperature criticality 

Consider a general perturbation H' which augments the 
Klein-point Hamiltonian to give 



Htot — H + H': 



(47) 



Eq. ( |23] | provides one specific example. We employ the no- 
tation e (Sec. IV) to denote the small energy scale of H' . In 
what follows, we prove that unless the ground states of Htot 
are ordered, in the sense that the different ground states are 
distinguishable by local order parameters, then a region of 
finite-temperature criticality must occur if the perturbation is 
not linear in an order parameter 

We begin by stating that the system displays critical cor- 
relations (i) when e = at any finite temperature and (ii) 
at T/e oo. We also note that (iii) any local order pa- 
rameter vanishes for any fixed e at sufficiently high temper- 
atures T > Tc{e). Finally we assume (iv) that at suffi- 
ciently high temperature the two-point correlator of an or- 
der parameter decreases monotonically in temperature at long 
enough distances. While (i), from Sec. IV, and (iii) are self- 
evident, the demonstration of statement (ii) represents the core 
of the proof: because critical correlations appear in the high- 
temperature limit, (iii) and (iv) imply that correlations at any 
finite but high enough temperature cannot be weaker than al- 
gebraic. Thus the system exhibits critical correlations at suffi- 
ciently high but finite temperature. 

The derivation of statement (ii) proceeds from the definition 
of the relevant order parameter By analogy with the classical 
dimer polarization operators presented in Sec. IIIC, we define 
the corresponding quantum operators. For a given plaquette A 
formed by the four sites (a, b, c, d) appearing in counterclock- 
wise order, the polarization is given by 

Pa = Pa-xCx + PA-.y^y + Pa-zCz 

= {Sa ■ Sd — Sb ■ Sc)ex + {Sa ■ Sb — Sc ■ Sd)ey 

+ {Sa-S,-Sb-Sd)e,'\. (48) 

The dimer-dimer correlation function between two crossed 
plaquettes A and B is then given, under the assumption that 



all states in the Klein manifold are linearly independent 14511 . 
and at finite temperatures T e, by the expectation value 

XAB = {Pa ■ PB)<L.T>a 



Tr 


'e-P"PA ■ 


Pb 


Tr 





^lim XAB = 



oc \rAB\ 



(49) 



with \rAB \ the separation between the centers of the plaque- 
ttes A and B, d the spatial dimensionality, Ng the number of 
states in the ground-state manifold (Sec. IIIA) and the sum 
taken over a complete set of orthonormal states {iV")}- At 
the Klein point, the correlation function xab is algebraic, as 
shown in Sec. IV. 

The essential property underlying the physics of the finite- 
temperature derivation is that the value of xab is the same for 
all temperatures when the magnitude e of the external per- 
turbation vanishes. This is a consequence of the fact that 
all states within Klein-point manifold carry equal weight, 
whence xab is given by the last two lines of Eq. ( |49l ) and 
the system is manifestly critical. As the Klein-point limit 
is approached for T/e oo, this average becomes equiva- 
lent to a sum over all states spanning the Klein-point Hilbert 
space and its value is precisely the same as for the T = 0+ 
case when e is taken to zero. Here the probability of any 
state \tp) as given by the finite-temperature density matrix is 
exp[—/3E^]/Z = 1/Ng. For T/e — > oo the correlation func- 
tion XAB then takes the form specified by Eq. ( |49] l, decaying 
algebraically to as stated in (ii) above. 

Although it is not relevant to the above derivation, we re- 
mark briefly that topological rules similar to those of Sec. V 
may be devised for the computation of the dimer-dimer corre- 
lation function. As a consequence of the two scalar products 
appearing in each term of xab, two permutation operations 
can be performed on the loop formed by two states (as op- 
posed to the single operation discussed in Sec. V). A similar 
set of topological rules for the computation of xab will be 
provided elsewhere. 

The use of statements (ii) and (iii) for arbitrary pertubations 
is assisted by considering the form of general order parame- 
ters. Let {Ipq)} denote the ground states of Htot and O the 
order parameter by which they are distinguished. In a system 
which displays local order. Or is an operator which depends 
on a finite number of fields in a sphere (or disk) of finite radius 
R which surrounds a given point r). A general algorithm for 
the construction of the order parameter by examining reduced 
density matrices is given in Ref. [68[|. 

At T = within any ground state, (OrOr') ^ |(Or)P = 
^ for separations |r — r'| ^ oo. Following Eq. ( |49l ). 



(Or 



E^,|<^.)(V'l</>l)Mi2(02|a|03)M34(04|^>e-'3^'* 



(50) 
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where denotes a set of pure dimer states and a 

complete set of orthonormal states spanning the Klein-point 
basis. Equation dSOl l makes use of the overcompleteness rela- 
tion 

Klein, (51) 

01,02 

with iKicin = 11^(1 ^ '^he unit operator in the Klein- 
point subspace [and the projection operator onto the sub- 
space of net spin ~ 2 (Sec. II)]; Eq. ( fsTl ) expresses the fact 
that the dimer states completely exhaust the Klein-point sub- 
space, as proven in Ref. |38]. As emlpoyed in ll45ll . all Klein 
model ground states can be expressed as a sum of projection 
operators on these ground states which decompose them into 
(Young tableaux) sectors which are well defined under the 
permutation of pair sites. The presence of a singlet dimer 
between two sites enforces antisymmetry amongst the perti- 
nent lattice sites. As evident in expressing the ground states 
in terms of the irreducible representations of the permutation 
group, the most general operations allowed within this basis 
amount to permutations. ll69ll As the most general permuta- 
tion can be written as a product of pairwise permutations and 
as Si ■ Sj = \{Pij ~ where is the operator permut- 
ing sites i and j, the most general permutation is a functional 
of scalar spin products. For an operator Or which is a func- 
tional of spin products of the form Sa ■ Sb, the correlation 
function (Or Or') is equivalent to a correlation function be- 
tween dimer products. From the analogy to a dipolar system 
(Sec. Ill), correlation functions containing more scalar spin 
products decay algebraically according to the strength of the 
Coulomb interaction between more complicated dipole con- 
figurations. Moments of higher order than dipolar are impor- 
tant for those dimers which are close to each other but far 
from all other dimers. All correlation functions thus decay 
algebraically with the dimer-dimer separation. The analysis 
based on Eq. ( |49] l may now be extended to more general corre- 
lation functions which involve the order-parameter fields Or. 

Returning to statement (iii), as in our earlier considerations 
for the polarization P this is based on the result that any 
finite-temperature expectation value (Or Or') will decrease 
monotonically with increasing temperature, and asserts that 
when the expectation value of the order parameter vanishes, 
(Or) = 0, the corresponding energy scale sets the transi- 
tion temperature [(iii)]- This quite universal property is 
reviewed briefly in Ref. [|72|]. We reiterate [(ii)] the fact that 
in the high-temperature limit, where correlations are weak- 
est, the correlation function (OrOr') retains a decay which is 
only algebraic in the separation |r — r'|. From these results 
it is clear that there is a finite-temperature region surrounding 
the "Klein line" (e = 0, T > 0), and specified explicitly by 
T > T'c(e), where critical correlations appear 

We conclude the general discussion by noting that for a 
perturbation which is linear in a general order parameter, 
H' = —eJ^A Qa, it is possible that the two-point correla- 
tion function takes a finite value, (QaQb) fn^ 7^ as 
I^abI 00, at any finite T. This type of perturbation may 
then yield correlations which are algebraic only in the limit 



T ^ 00 iItsIi . In the next subsection we consider a partic- 
ular perturbation whose effect can be approximated by the 
anisotropic six-vertex model to show that the dimer correla- 
tions may also remain critical up to a threshold value of the 
perturbation. 



B. Anisotropic perturbation 

In this subsection we use a specific anisotropic perturba- 
tion to investigate the nature of the phase diagram. We will 
illustrate explicitly that the phase diagram contains a finite- 
temperature critical region, exemplifying the general concepts 
of the preceding subsection, and then derive its approximate 
phase boundaries. 

We consider a perturbation of the form 

Hd = Jd ^ Si ■ Sj, (52) 

where {{ij)) denotes all diagonal pairs on the same crossed 
plaquette. This type of interaction emulates the weakening 
(for Jd < 0) of exchange interactions between neighboring 
diagonal sites on a checkerboard plaquette relative to those of 
horizontal or vertical spin pairs, or a lifting (to tetragonal) of 
the cubic symmetry in a pyrochlore tetrahedron. Focusing for 
specificity on the checkerboard geometry, for Htot = H + Hd 
it is clear that in the regime Jd,T <^ J, the effect of Hd 
is to favor certain dimer, or equivalently six-vertex, config- 
urations over others, and in particular that for Jd > the 
diagonal contribution, {tl.'a\Hd\il.'a), favors singlet formation 
on the diagonal bonds of all plaquettes. Because off-diagonal 
brackets are significantly smaller than diagonal contributions, 
these may be neglected. The favored ground states in this case 
are precisely the two antiferroelectric states. Similarly, for 
Jd < the single dimer on each crossed plaquette lies pref- 
erentially on the vertical or horizontal bonds. Because only 
the diagonal contributions {il-'a\Hd\tpa} are considered, the 
ground-state degeneracy for Jd > remains exponentially 
large: the system supports that subset of all allowed ground 
states which corresponds to the ground-state manifold of the 
Klein model on the square lattice, the size of which scales 
with 2^+^, where N is the number of lattice sites We 
remind the reader that the states in this restricted set of dimer 
coverings, in which precisely one dimer is present on every 
square-lattice plaquette, are strictly orthogonal to each other 
in the thermodynamic limit. The general case has a very much 
larger manifold of ground states, the classification of which 
proceeds by the line representation presented in Sec. III. 



1. Critical region 

The illustration that the phase diagram in the presence of 
this perturbation contains a finite-temperature critical regime 
is based on the exact limits examined in Subsec. (I VI At and 
constitutes an explicit implementation of this discussion. We 
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begin by noting, as before, that when = the system re- 
mains at the Klein point for all temperatures, displaying crit- 
ical correlations, and the analysis of Sec. IV is applicable. 
The susceptibilities x'ab = {PA-aPB-.a) are a direct probe 
of dimer correlations: at = the correlation function xab 
is critical for all T > 0. 

With reference to the considerations below Eq. (|52] |. for 
Jd > 0, Xab increases in magnitude relative to its Klein- 
point value for diagonal dimers only. Thus, ( i) the decay of x 
with distance is no faster than algebraic. However, no VBC 
order is possible at high temperatures T ^ Jd, and thus ( ii) 
X — > as the plaquette separation |r| oo. Conditions 
( i) and ( ii) taken together imply that for T ^ Jd the system 
exhibits algebraic correlations. Similar considerations apply 
when Jd < 0, where [|xabI + I Xab I] increases relative to its 
Klein-point value, corresponding to correlations of horizon- 
tal or vertical dimers. We thus conclude that for T ^ \Jd \ the 
system is critical. These statements are completely general, 
and provide one explanation for why it is that the six-vertex 
model always has critical, and only critical, correlations at suf- 
ficiently high temperatures. ll74ll 



2. Finite-temperature criticality by approximate mapping to an 
anisotropic six— vertex model 

We now derive the approximate forms of the phase bound- 
aries of the critical region. This calculation is not exact for 
the reasons stated at the start of this section, and is included 
to provide a qualitative indication of the nature of the phase 
diagram. We begin with an explicit statement of the approxi- 
mations employed in this derivation, (i) In contrast to Sec. V, 
for this finite-temperature derivation the non-orthogonality 
of the different dimer states is neglected. At the Klein point 
{Jd = 0) all dimer configurations carry equal weight and their 
orthogonality is irrelevant, a fact exploited in Sec. V. All of 
the properties of the system are a consequence of the entropic 
considerations presented in Sec. IV. However, away from the 
Klein point the configurations no longer have equal proba- 
bility (in the example of Eq. ([52]), states with diagonal dimer 
configurations are favored when Jd > 0, and suppressed when 
Jd < 0). (ii) Off-diagonal contributions in the full Hamilto- 
nian, i.e. terms of the form {tpa\H\tpi,) with a 7^ 6, are ne- 
glected. 

When mapped to the six-vertex representation, provided 
that the conditions Jd,T ^ J are satisfied, the anisotropic 
perturbation of Eq. ( |52b is reflected in a change in the relative 
weights of the six-vertex configurations (which at low tem- 
peratures are the only states with non-vanishing weights). At 
T = 0, for Jd > these nonuniform weights are represented 
by vanishing energies for four of the vertex types while the 
other two contribute an equal and finite energy, and conversely 
for Jd < 0. 

With reference to Fig. [3] we choose the energies of the 
six-vertex configurations displayed to be ei=i,2,3,4 = and 
£5 = £6 = ^l^d in the presence of the perturbation. The cor- 
relations and thermodynamics of this model are best analyzed 
through the fugacities a ~ cxp[— /3ei] = exp[— /3e2], b — 



T' 






y AFE-VBC 
f- > 









J. 



FIG. 17: Approximate phase diagram for Klein-point Hamiltonian 
with diagonal exchange perturbation Jd- Neglecting both non- 
orthogonality and off-diagonal contributions, an antiferroelectric va- 
lence bond crystal (AFE-VBC) ordering is stabilized below T = 
3Jd/(81n2) for Jd > 0. The system is critical at all points within 
the region shaded in light grey. In the absence of off-diagonal pro- 
cesses, it remains critical to zero temperature for Jd < 0; the qualita- 
tive effect of these processes is to remove the low-temperature criti- 
cal behavior by stabilizing a VBC ordering in the region shaded dark 
grey. While the phase boundaries shown above are approximate, the 
existence of a finite-temperature critical region is not (Sec. VIA). 



cxp[-/3e3] = exp[-/3e4], c = exp[-/3e5] = exp[-/?e6], 
where (3 is the inverse temperature. The fugacities determine 
a parameter A whose value determines in turn the phase in 
which the six-vertex model lies. 



62 



2ab 



(53) 



Following Lieb's original solution of the six-vertex model 
Ii42i i43il , all models with |A| < 1 are in the disordered phase. 
Such small-| A| systems, including the ice model relevant for 
the isotropic checkerboard lattice where a = 6 = c = 1, are in 
the disordered phase and exhibit power-law coiTelations. By 
contrast, for A > 1 the system adopts an ordered configura- 
tion where all arrows have a definite chirality (such as up and 
to the right or down and left). In the system with perturbation 
Jd, a = b = 1 while c = cxp[|/3Jrf], and one finds that if 



T > 



Sin 2 



Jd, 



(54) 



then |A| < 1, indicating that the system is in its disordered 
phase. In this case we find power-law correlations among 
dimer pairs and the system is trivially critical for all non-zero 
temperatures when Jd < 0, or critical only at sufficiently high 
temperatures [whose lower bound is given by Eq. ( l54l il when- 
ever Jd > 0. When T = 0, one finds that for all Jd < 
the system lies precisely on the boundary between the ordered 
and the critical phase (A = 1). Equation (l54l i defines a finite 
region within the {Jd,T) plane (Fig.[T7]i in which the system 
is critical. 

The above results are applicable under the conditions that 
all relevant states lie within the Klein basis manifold, i.e. T <C 
J, that the effects of non-orthogonality are neglected, and that 
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only the diagonal corrections induced by Hd are considered. 
Although these are small, off-diagonal contributions of the 
type {ilJa\Hd\^b), analogous to those computed in Sec. V, are 
relevant in removing the "marginally critical" T = Q corre- 
lations obtained for Jd < 0, and stabilize once again a par- 
ticular valence-bond ordering pattern of the types considered 
in detail in Sec. V. In this situation the marginally critical line 
may be moved to finite temperatures, separating the low-T 
ordered phase from a finite-T critical regime. In summary, 
it is possible to demonstrate rigorously that specific types of 
perturbation around the Klein point lead to an exotic, exactly 
critical regime. Quantitative calculations of the properites of 
this regime are not exact, but approximations of the kind ap- 
plied here may be used for the qualitative determination of 
quantities such as the boundaries of this phase. 




Be B 

FIG. 18: (Color online.) Low-temperature phase diagram for the 
Klein Hamiltonian in the presence of a Zeeman term (applied mag- 
netic field B) for the checkerboard and pyrochlore spin systems. 



VII. DISCUSSION AND CONCLUSIONS 

The pyrochlore (A2B2O7) and spinel (AB2O4) structures 
are very common among magnetic oxides, occurring for a 
wide variety of transition and rare-earth metals as both the 
A and B ions and thereby offering a range of valence states 
and magnetic moments 5 in a pyrochlore lattice of interacting 
spins. Because the orbital contribution to the magnetic mo- 
ment of transition metal ions in the 3d series is quenched by 
the crystal-field splitting, the spin-orbit coupling has only a 
perturbative effect and the relevant physical spin models are 
effectively SU(2)-invariant. While the link to a physical py- 
rochlore or spinel system satisfying all the necessary criteria 
remains to be found (Sec. II), we have shown here that for 
these structures a class of SU(2)-invariant 5 = 1/2 models 
of this type, whose local interactions emerge directly from 
the simplest Hubbard Hamiltonian, gives rise to spinon ex- 
citations which propagate in the full lattice. 

At the root of the exotic behavior exhibited by this type of 
system is the massive (extensive) degeneracy of the ground- 
state manifold in the vicinity of the Klein critical point. We 
have gained physical insight into the multiple essential as- 
pects of this problem by exploiting exact mappings between 
the spin system under consideration and other physical sys- 
tems [75]. In particular, because the ground-state sector of 
the 5=1/2 Klein models on the pyrochlore lattice obeys ice 
rules of the type discussed in Sec. Ill, it is possible to exploit 
analogies with six- vertex models p3| . string gases and U(l) 
gauge magnets to make a number of powerful qualitative and 
quantitative statements. The local constraint of zero diver- 
gence imposed by the ice rules implies that the low-energy 
sector is described by an effective U(l) gauge theory in the 
neighborhood of the Klein point. This theory is, however, 
not necessarily the minimal U(l) gauge magnet, whose terms 
involve only the smallest loops, which is usually invoked in 
these cases ll20ll . The variational approach which we intro- 
duce and employ in Sec. V indicates that in some cases, as 
a result of the non-orthogonality of the singlet dimer cover- 
ings, processes involving larger loops are relevant for realistic 
physical perturbations. 

Although the spinons are deconfined excitations only at a 



single point (the Klein point) of the quantum phase diagram 
(T = 0), this deconfinement exists over a finite region of the 
phase diagram due to thermal fluctuations which give rise to 
a dilute plasma phase of spinons. The spinon-spinon corre- 
lations in this region are dominated by the T — > 0+ Klein 
critical point, whose critical fluctuations are of entropic ori- 
gin. The primary characteristic of the Klein point is the sta- 
bilization of a low-energy manifold of states satisfying a lo- 
cal constraint of zero divergence, or equivalently the ice rules 
ifTolfTlll . and the removal of any low energy scale. As a con- 
sequence of both conditions, the system becomes critical due 
to entropic fluctuations induced even at infinitesimal temper- 
atures. 

These features are directly relevant for the potential ob- 
servation of spinons in two- or three-dimensional systems. 
In the case of the pyrochlore lattice models considered here, 
the deconfined phase is expected to appear near a transition 
between two valence-bond crystals. Similar considerations 
were applied to the S* = 1/2 model studied in Ref. the 
difference in that case being that the spinons propagate along 
one-dimensional paths because the dimensional reduction at 
the corresponding Klein point is incomplete. 

It is also interesting to note that an exact deconfined quan- 
tum critical point (QCP) jst] can be induced by adding a Zee- 
man term to the Klein Hamiltonian of Eq. (O. The QCP 
is induced at the critical magnetic field Be which closes the 
singlet-triplet spin gap. The excited two-spinon deconfined 
state becomes the ground state at Be. For higher values of 
the field {B > Be), the concentration of spinons becomes fi- 
nite and the degeneracy between the underlying singlet dimer 
coverings is lifted. Consequently, spinon-spinon confinement 
is restored and the resulting triplet pairs are expected to con- 
dense, giving rise to an XY-type antiferromagnet (XY-AF) in 
the plane perpendicular to the applied field. 

In the presence of a magnetic field, the system at the Klein 
point displays a true deconfined QCP between a line of clas- 
sical critical points in the region B < Be and a magnetically 
ordered region at B > Be with two coexisting order parame- 
ters. A schematic phase diagram illustrating the finite extent 
of the deconfined regime is shown in Fig. [18] Very different 
and much more conventional behavior is obtained away from 
the Klein point, where the ground-state degeneracy is lifted. 
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Here the region B < Be has a VBC ground state with gapped 
magnon excitations, separated by a conventional QCP from an 
ordered phase with spin-wave excitations. 

The ideas and concepts connected with deconfinement at 
quantum critical points have formed an extensive recent liter- 
ature best summarized in Ref. Jst]. As noted in Sec. I, many 
of these studies depart from effective U(l) gauge-theoretical 
treatments rather than microscopic models, and it is not yet 
clear that specific systems exist which realize the desired de- 
confinement properties. The Klein point in certain pyrochlore 
lattice models was considered briefly in Ref. ifsoll . although 
these authors did not dwell on either the detailed physical 
properties of such a point or on their origin. We have found 
that the mapping used in this study to a quantum dimer model 
with only RK interactions is not in fact justified for a Heisen- 
berg antiferromagnet in a pyrochlore geometry. The sce- 
nario of the "constrained entropic critical point" introduced 
recently in Ref. j?^ contains certain parallels to the physics 
of the Klein point with regard to high-temperature constraints. 
However, this is said to be a "top-down" construction based 
on different energy scales, some of which are present in a 
number of the effective models considered, and does not con- 
tain specific properties emerging from a microscopic Hamil- 
tonian for all temperature regimes. 

We conclude this discussion with a speculation. The ex- 
ponential ground-state degeneracy [Eqs. ( 1121141 )1 and result- 
ing extensive configurational entropy which we have found 
suggest, but certainly do not mandate, that glassy spin dy- 
namics may occur naturally in pyrochlores. We stress that 
these systems are uniform and disorder-free but frustrated, 
whence this suggestion would reinforce the widely held be- 
lief that frustration rather than disorder is the fundamental re- 
quirement for the dynamical properties of structural glasses. 
The basic premise of the Vogel-Fulcher (VF) form of glassy 
dynamics is that the relaxation times appear to diverge at a 
temperature which correlates well with the intercept (Tq) of 
the extrapolated entropy of the super-cooled liquid to that of 
the solid. According to a prevalent line of reasoning iFttIi . 
as the system is cooled the configurational entropy first be- 
comes extensive at Ta, the onset temperature of multiple local 
free-energy minima; at a lower temperature, Tq, these min- 
ima become stable and the configurational entropy vanishes. 
If an exponentially large number of metastable states is found 
for Ta > T > Tq, one may invoke the analysis of entropic 
droplets IttIi to obtain the characteristic free-energy barrier 
height 

AE {TScy^. (55) 
Linearizing the extensive configurational entropy, 

5e(r ^ r+) ~ V{T/To - 1), (56) 
leads to VF dynamics iIttIi with relaxation times 

T^exp[DTo/iT-n)]. (57) 

Here Tn a temperature scale close to the Kauzmann temper- 
ature IITSIl at which an "ideal glass transition" would occur. 



and where the extrapolated entropy of the assumed liquid un- 
dergoes a crisis. We stress that this very general derivation 
of a systematic definition for glassy behavior, by consider- 
ation of free-energy barriers iIttIi . does not require that the 
different ground states be linked only by infinite-length pro- 
cesses. (However, if only the spatially longest processes were 
operational this would indeed lead to very robust, exponen- 
tially slow quantum dynamics, as first noted by Chamon and 
coworkers ll79ll .) For the pyrochlore system we have shown 
that there exist finite-length tunneling processes linking dif- 
ferent ground states within the highly degenerate low-energy 
sector, and thus suggesting an obvious candidate system for 
this type of glass. We comment that the physical processes 
indicated here, namely of glassy spin dynamics arising from 
frustration in a periodic lattice, are similar to early ideas of the 
"topological spin glass" of Ref. I8O1I . Empirically, glassy dy- 
namics are strongly indicated in a number of frustrated mag- 
netic systems, of which we quote only a selective list: evi- 
dence of phases with certain glassy characteristics has been re- 
ported in some pyrochlores, including Y2M02O7 1181 118211 . in 
the stacked kagome layer system SrCr8.gGa3.4O19 Il82ll . and 
in the kagome bilayer compound Ba2Sn2ZnCr6.8Ga3 2O22 
llssll . Further indications for glassy behavior are found in 
some quasi-triangular antiferromagnets, including NiGa2S4 

a. 

In conclusion, we have performed a detailed analysis of 
a physically motivated quantum spin model with only near- 
neighbor interactions on the two- and three-dimensional py- 
rochlore lattices. This reveals a wealth of exotic behavior 
which can be traced to the extensive degeneracy of the sys- 
tem at the Klein point and thus to a complete dimensional re- 
duction. The complete absence of quantum fluctuations in the 
ground-state manifold at the Klein point leads to a new type 
of classical critical point, with exactly known critical behav- 
ior, which we stress emerges from a highly frustrated, micro- 
scopic quantum spin model. 

The critical correlations at the Klein point are driven by en- 
tropic (thermal) fluctuations, leading to an effective plasma 
phase with Coulombic interactions. These are not confin- 
ing, whence the elementary excitations are deconfined spinous 
which propagate freely in all directions (a further mani- 
festation of complete dimensional reduction). This high- 
dimensional fractionalization would also be manifest as a 
spin-charge separation in the dilute limit of added holes if 
the carrier hopping were much smaller than the magnetic en- 
ergy scale J. Physically relevant perturbations away from 
the Klein point lead in general to confined phases of static 
valence-bond order with spin gaps to 5 = 1 excitations. 
However, at finite temperatures in the vicinity of the Klein 
point, the classical criticality is dominant and deconfinement 
persists over a finite region of the phase diagram. This type 
of deconfined behavior goes well beyond the simple, thermal 
decoupling of a system (to obtain quasi-one-dimensional be- 
havior) in that the classical critical point is driven to T = 0+: 
the microscopic origin of deconfinement lies in the proximity 
of the system to a Klein point, and spinon propagation remains 
d-dimensional. 

We have obtained our conclusions, and as a result a rather 
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complete picture of a distinctive paradigm for classical and 
quantum criticality, from a number of rigorous techniques. 
The dimer coverings of the checkerboard and pyrochlore lat- 
tices may be treated, through their connection with the ice 
rules, by both the six-vertex mapping and a line represen- 
tation, which exploit the underlying topological order of the 
physical system to classify the states of the manifold. We have 
developed a loop representation of the physical processes in 
the dimer basis which has very general applicability to calcu- 
lations involving the non-orthogonal dimer states; this gives a 
clear and intuitive picture of all contributing local processes, 
reflected in the sizes of corresponding loops, and a straightfor- 
ward but rigorous set of rules systematizing their computation. 
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Note added in proof 



Two months after the initial appearance of our work f85\, 
Pollmann and coworkers reported, in a series of preprints [18611 . 
on a spinless Fermi model which possesses several of the qual- 
itative features which we found here for spin systems on the 
pyrochlore and checkerboard lattices. 
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